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Abstract — In this paper, a family of low-density parity-check 
(LDPC) degree distributions, whose decoding threshold on the 
binary erasure channel (BEC) admits a simple closed form, 
is presented. These degree distributions are a subset of the 
check regular distributions (i.e. all the check nodes have the 
same degree), and are referred to as p-positive distributions. It 
is given proof that the threshold for a p-positive distribution 
is simply expressed by [A'(0)p'(l)]^. Besides this closed form 
threshold expression, the p-positive distributions exhibit three 
additional properties. First, for given code rate, check degree and 
maximum variable degree, they are in some cases characterized 
by a threshold which is extremely close to that of the best known 
check regular distributions, under the same set of constraints. 
Second, the threshold optimization problem within the p-positive 
class can be solved in some cases with analytic methods, without 
using any numerical optimization tool. Third, these distributions 
can achieve the BEC capacity. The last property is shown by 
proving that the well-known binomial degree distributions belong 
to the p-positive family. 

I. Introduction 

The unavailability of a closed form expression for the 
decoding threshold of low-density parity-check (LDPC) codes 
still represents an open problem. So far, no general closed 
form threshold expression is known for any transmission 
channel. Some results in this sense have been developed for 
the binary erasure channel (BEC), for which some analytic 
threshold expressions have been proposed. In [1], an analytic 
threshold expression has been presented for regular LDPC 
ensembles. A more general analytic expression, valid for check 
regular ensembles, has been proposed in [2]. To the best of 
the authors' knowledge, the most general analytic expression 
currently available is that one developed in [3], which can 
be applied to fully irregular ensembles. For regular or check 
regular ensembles, the formula proposed in [3] coincides with 
that one from [2]. The problem with these expressions is that 
none of them can be really considered a closed form one 
(except for the case of regular LDPC ensembles with degree-2 
variable nodes presented in [1]). In fact, in all these formulas, 
the threshold is a function of some parameter which in general 
does not admit a closed form expression. This parameter, 
which depends on the degree distribution, can be a root of a 
real polynomial, a fixed point of a real function, or the abscissa 
of the tangent point between two EXIT curves [4]. 

An exception is represented by the degree distributions for 
which the first occurrence of a tangency point between the 



EXIT curves in the EXIT chart appears for a value of the a pri- 
ori mutual information equal to 1 (this is sometimes referred 
to as derivative matching condition). Denoting the decoding 
threshold by 5* , these distributions achieve with equality the 
stability condition 5* < [A'(0)p'(l)]"^ [5] (which holds for 
any distribution). This paper presents a family of check regular 
LDPC distributions, called p-positive distributions, which are 
a subset of the distributions achieving the derivative matching 
condition. The starting point for obtaining this family is the 
afore mentioned analytic formula developed in [2]. 

The name "p-positive" is chosen because these distributions 
are defined as those ones for which a certain polynomial 
p (•), whose coefficients depend on the edge-oriented variable 
distribution X{x) and on the degree dc of the check nodes, is 
non-negative between and 1 . The theory of polynomials can 
be applied to obtain conditions for a check regular distribution 
to be p-positive. A useful theorem in this sense is the Fourier- 
Budan theorem [6, p. 27], which states an upper bound to the 
number of real roots in a given interval, for a polynomial with 
real coefficients. The application of the Fourier-Budan theorem 
to the polynomial p (•) leads to a set of inequalities involving 
the coefficients A^'s and dc, which represent a necessary 
condition for a distribution to be p-positive. 

Besides the closed form threshold expression, the p-positive 
distributions are shown to exhibit some additional good prop- 
erties. First, the threshold of the optimal p-positive distribution 
under a set of constraints, represented by the given code rate, 
check degree and maximum variable degree, is in some cases 
extremely close to that of the best known check regular degree 
distributions. Second, within the p-positive class, it is in some 
cases possible to optimize the degree distribution with analytic 
methods only, i.e. without using numerical optimization tools. 
In fact, despite being only necessary, the condition obtained 
from the Fourier-Budan theorem is shown to be sufficient, in 
some cases, to find the p-positive distribution with optimum 
threshold under the given set of constraints. Third, it is proved 
that capacity achieving sequences exist within the p-positive 
class. More specifically, it is recognized that any binomial 
degree distribution [7] is p-positive. 

The paper is organized as follows. Section recalls the 
threshold expression for check regular ensembles from [2]. 
Section HU] exploits this formula for introducing the new class 
of distributions, and develops a necessary condition for a 



check regular distribution to belong to this class. Section 
HVl is devoted to the threshold investigation of the p-positive 
distributions, with both differential evolution [8] and analytic 
methods. In Section |V] the existence of capacity achieving 
sequences within the p-positive family is proved. Section [Vl] 
concludes the paper. 

II. Threshold for Check Regular Distributions 

Let C"(A, p) be the ensemble of all the length-n LDPC 
codes with edge-oriented degree distribution (A, p) [5]. Let 6 
be the BEC erasure probability. Then, the (bit oriented) thresh- 
old 6* for the ensemble C^{X,p) on the BEC is defined as 
the maximum 6 for which the residual bit erasure probability 
can be made arbitrarily small by increasing the number of 
decoding iterations, in the limit where the codeword length n 
tends to infinity. 

Density evolution for the BEC [5], can be expressed as 
follows. If Xi is the probability that a message from a variable 
node to a check node during the ^-th decoding iteration is an 
erasure message and the bipartite graph is cycle-free (infinite 
codeword length), then x^+i = 6 \{1 — p{l — xi)). Hence, 
the threshold 5* is the maximal 5 for which in the 

limit where ^ ^ oo. The first appearance of a non-zero fixed 
point for xt is a necessary and sufficient condition for the 
corresponding 5 to be the threshold. 

From this observation, 5* can be equivalently expressed as 
the maximal 5 for which 5 A(l — p(l — x)) < x V x E {0,6] 
[9], or the maximal S for which: 



p{l-S\{x)) >l-x Vx G (0,1] 



(1) 



For check regular codes, with variable distribution X{x) = 



and check distribution p{x) 



(dc > 3), 



the inequality ^ becomes: 

(1 - S\{x) f--'^ > 1 



X Vx e (0,1]. 



(2) 



Let us define f{x,6) = (1 — dX{x))'^'=^^ and suppose at 
first A2 — 0. For any given x G (0, 1], f{x,S) is continuous 
and monotonically decreasing with respect to S, varying from 
/(x,0) = 1 to fix, I) = (1 - A(.t))'^--i > 0. Moreover, 
it is everywhere derivable for x in (0,1), and the following 
relationships hold: 

df/dx{Q,S) =0 V(5 G (0,1) 
f{0,S)= 1 V<5g(0,1) 
f{l,5)= {l-5f'-'>0 
f{x,0) = 1 VxG (0,1). 

It follows from these properties that the maximal value of 6 
for which (|2]l holds is the minimum value of 6 for which the 
graph of f{x, S) (considered as a function of x with 6 playing 
the role of parameter) is tangent to the graph of g{x) = 1 — x, 
for some a; = 7. The tangency condition in x = 7 is: 



S{d, - 1)A'(7)(1 



1-7 

- sx{j)r^- 



(3) 



The relationships ^ are not able to unambiguously de- 
termine the threshold. In fact, (O can admit several discrete 
solutions (5,7), with 7 G (0,1) and 5 G (0,1). From a 
geometrical point of view, several discrete values of (5 G (0, 1) 
can exist for which the graph of f{x, 5) (for fixed 5) is tangent 
to the graph of (7(x) = 1 — x in some x — j. The threshold 
S* for the check regular distribution is the minimum among 
these discrete values of S. 

The second equation of (O can be written as 



(1 - 5Xm 



[5{d, - 1)A'(7)]- 



that can be substituted in the first equation in order to obtain 



<5=[A(7) + (4-l)(l-7)A'(7)]- 



(4) 



In the following, h{-) will denote the following key function: 

h{x) = [X{x) + {d, - 1)(1 - x)X'{x)]-\ (5) 
By substituting ^ into the first equation of ([3]), we obtain 

A(7) + (4-l)(l-7)A'(7) 



(de-l)(l-7)A'(7) 
and developing this expression leads to: 

A(7) 



(I-7) ''c-l. 



7 = 1 



(4 - 1)A'(7) 



de-1 

1-7)'^". (6) 



Summarizing, for a given check regular ensemble with A2 = 
0, the threshold on the BEC is equal to the minimum ^1(7) 
under the constraint that 7 is one of the (usually) several fixed 
points in (0, 1) of 



(x) = 1 - 



X{x) 



[dc - l)A'(x) 



dc-1 

l-x]'^-\ (7) 



Note that 7 = is always a fixed point for ((>{■). However, for 
A2 = the threshold is never achieved in correspondence of 
7 = 0, because in this case lima:^o+ h{x) = +00. 

The only differences when removing the hypothesis A2 = 
ai-e that df/dx{0,S) < V(5 G (0, 1), and that h{0) is finite 
and equal to [A2(4 - 1)]"^ = [X'{0)p'{l)]-^. Hence, the 
threshold could be achieved in correspondence of the fixed 
point 7 = 0. Thus, the threshold for a check regular ensemble 
with A2 > is equal to the minimum h{'j), with 7 fixed point 
in [0,1) of 0(.). 

In the special case of regular LDPC codes, with variable 



distribution X{x) 



dv > 3, and check distribution 



p{x) — x'''" ^, dc >'i, the threshold is given by 

[(4-i)(d. -1)]-^ 



(8) 



where 7 is the unique fixed point in (0, 1) of 

i/'(x) = 1 - (1 - cx)^ . (9) 
with c = [{dc - l)(d. - 1) - l]/[{dc - l)(rf„ - 1)] < 1. 



III. The New Class of Degree Distributions 

In this section, the class of p-positive degree distributions 
is introduced. 

Definition 1: A degree distribution with code rate R and 
threshold 5* = (1 — e)(l — i?) is called (1 — e) capacity 
achieving of rate R. 

Consider a check regular distribution for which h{x) > 
h{0) Vx G (0,1], where h{-) is defined in (|5]l. This can be 
equivalently written as < [/i(0)]^^, i.e. 

L L 

^A,x'^i + (4-l)(l-2:)^(i-l)A,x'-2 < A2(dc-1), 

i=2 i=2 

where L denotes the maximum variable degree. After some 
algebraic manipulation, the previous inequality can be put in 
the form p{x) > 0, where p{-) is the real polynomial 

L-l 

p (x) = UJL Al x^^^ + ^ [ujiXi - {uji+i + 1) Aj+i] a;*"^, 

i=2 

(10) 

where uji = (dc — 1)(* — 1) ^ 1- The polynomial p(-) will 
be expressed in the form p{x) = ^f^Q PiX^- The condition 
h{x) > h{0) for all x G (0, 1] is equivalent to the condition 
that the real polynomial p{-) is positive or null in (0,1]. 
The class of degree distribution pairs studied in this paper 
is introduced next. 

Definition 2: A p-positive distribution is any check regular 
degree distribution with A2 > 0, such that p{x) > V a; G 
(0,1]. 

The following theorem individuates a simple closed form 
for the threshold of the p-positive distributions, and states the 
necessary and sufficient condition for any p-positive distribu- 
tion to be (1 — e) capacity achieving of rate R. 

Theorem 1: The threshold of any p-positive degree distri- 
bution is equal to [A'(0)/9'(l)]~^. The degree distribution is 
(1 — e) capacity achieving of rate R if and only if 

A2 = [(l-e)(l-i?)(4-ir^ (11) 

Proof: For a p-positive distribution, h{x) > h{0) for all 
X G (0, 1]. From Section it is known that the threshold for 
a check regular ensemble with A2 > is the minimum among 
discrete values /i(7), with 7 fixed point in [0, 1) of </>(•) defined 
in (|7]i. Moreover, 7 = is always a fixed point of (p{-) under 
the condition A2 > 0. If h{x) > h{0) for all x G (0, 1], 
then the minimum is always achieved in correspondence of 
the fixed point 7 = 0, independently of the number and the 
positions of all the other fixed points of ((>{■). Hence, S* = 
h{0) = [A'(0)p'(l)]-i = [A2(4-l)]"^Ifand only if equality 
dlDholds, itis J* = (l-e)(l-i?). □ 

For a given e > 0, the search for p-positive (1 — e) capacity 
achieving of rate R distributions, with check degree dc and 
maximum variable degree L, can be performed by letting A2 
be expressed by (fTTT l. and looking for A;, i = 3, . . . , L, such 



that p{x) > for all x G (0, 1]. In general, for some R, dc 
and L, this problem admits solutions for e > Copt > 0, with 
Copt depending on R, dc and L. The value Copt is associated 
to the optimal p-positive distribution, under the imposed set 
of constraints (i.e. code rate, check nodes degree and active 
variable degrees). 

It is readily shown that p(l) — {dc ~ 1)A2 — 1, which is 
positive (from (fTTT l'). Then, the condition that p (x) > for x 
between and 1, is equivalent to the condition that p( ) has 
only roots with even multiplicity between and 1. Several 
well known properties of the real roots of polynomials can 
be then exploited, in order to develop conditions for a check 
regular distribution to be p-positive. In particular the following 
theorem, known as the Fourier-Budan theorem, permits to 
obtain a simple necessary (but not sufficient) condition. It 
provides an upper bound to the number of zeroes of a real 
polynomial between two values a and b, with a < b. 

Theorem 2 (Fourier-Budan Theorem [6, p. 21]): Let p(-) 
be a real polynomial of degree q, and let a and b be two real 
values such that a < b, and p{a)-p (b) ^ 0. Then, the number 
of real roots of p (•) between a and b (each one counted with 
its multiplicity) is not greater than A ~ B, where A and B 
are, respectively, the number of sign changes in sequences: 

p{a),p'{a),p"{a),...,p^'^\a) 
p(6),p'(6),p"(6), ...,p(')(6). 

Moreover, if the number of real roots of p (■) between a and h 
is smaller than A~ B, then it differs from A — B by an even 
number. 

This theorem can be directly applied to the polynomial p (•) 
defined by (fTOl i. where a = and 6=1, leading to the 
following corollary. 

Corollary 1: The number of real roots between and 1 
of p (x) = ^^^=0 Pi defined by ( fTol i is not greater than the 
number A of sign changes in the sequence po, p 1, . . . , p l-2- 
If the number of roots is smaller than A, then it differs from 
A by an even number 

Proof: The i-th derivative of p() in and 1 is 

p'^\<^)=^\p^ 

From the structure of the coefficients of p(-) it is not 
difficult to recognize that all the values in the sequence 
p (1), p'(l), • . . , p*'^~^-'(l) are positive for any L and dc > 3. 
Hence, 5 = and the number of roots of p(-) between 
and 1 is not greater than the number A of sign changes in 
the sequence po, l!pi, ■ ■ ■ ,{L ~ 2)!pi_2- This is equal to 
the number of sign changes inpo,pi,...,PL-2-By directly 
applying Theorem|2] we obtain the statement. We also observe 
that p i_2 — [{L — l)dc — L] Al is always positive. □ 

Recall that, for the p-positive distributions, p(-) has only 
roots, between and 1, with even multiplicity. Then, the 



following necessary condition for a check regular distribution 
to be a p-positive distribution is obtained from Corollary [T] 

Corollary 2: Necessary condition for a check regular dis- 
tribution to be a p-positive distribution is that the number of 
sign changes in the sequence po, P i, ■ • ■ j P l-2 is even. 

In the next section, some examples of threshold optimizations 
for p-positive distributions are provided. These results reveal 
that, for some values of the code rate R, check degree dc 
and maximum variable degree L, the p-positive distributions 
are characterized by very good thresholds. In some cases, 
the threshold is very close to that one of the best known 
distributions, under the same set of constraints. 

IV. Optimization of p-positive Degree 
Distributions 

A. Optimization with Constrained Differential Evolution 

The threshold optimization problem for the p-positive dis- 
tributions can be in principle performed numerically, using 
a constrained version of the differential evolution (DE) algo- 
rithm. Specifically, a number A^^i + l of values Xi = {l/ND)-i, 
for i = 0, . . . , Nd, are first chosen, for sufficiently large Njj, 
and then the DE optimization tool is run for given code rate, 
given check degree and given maximum variable degree, and 
with the additional set of constraints p [xi) > for all i. 

We performed this constrained DE optimization for R = 
1/2, L = 20 and 4 = {6,7,8}. Under the same set of 
constraints, we also performed the DE optimization, without 
imposing the p-positive bound. The results of this search are 
shown in Table H] for dc — 6 and dc — 7. In both cases, the best 
p-positive distribution exhibits a threshold that is only slightly 
worse than that one of the best check regular distribution ob- 
tained removing the p-positive constraint, with the advantage 
represented by the closed form threshold expression. 

This conclusion is no longer valid for dc = 8. In fact, for 
this check node degree, the best found p-positive threshold was 
S* = 0.469592, while the best found check regular threshold 
was S* = 0.491988. This means that the p-positive constraint 
is not "compatible" with all the possible sets of bounds on 
the code rate, check degree and maximum variable degree. 
When this compatibility holds, the p-positive distributions 
exhibit very good thresholds. In the other cases, they exhibit 
a threshold loss with respect to the best known distributions. 

B. Optimization Based on the Fourier-Budan Theorem 

In Section IIV-AI the p-positive distribution optimization 
has been performed according to a constrained version of 
the DE algorithm. Some examples of p-positive threshold 
optimization, based on a totally different technique, are pre- 
sented next. More specifically, it is shown that the necessary 
condition, developed in the previous section (Corollary |2]), 
can be exploited in order to perform the optimization process 
without using numerical tools. This is possible for particular 
structures of the variable nodes degree distribution. 

As a case study, let us consider check regular distributions 
with variable distribution in the form \{x) = A2 x + A3 x^ + 



TABLE I 

Optimal p-positive and not p-positive distributions and 

THRESHOLDS FOR i? = 1/2, L = 20 AND dc = {6, 7}. 





dc 


= 6 


dc 


= 7 




p-pos. 


not p-pos. 


p-pos. 


not p-pos. 


A2 


0.415884 


0.415273 


0.339162 


0.338843 


A3 


0.165968 


0.160268 


0.138401 


0.140058 


A4 


0.095028 


0.142202 


0.104711 


0.104198 


As 


0.106071 




0.033138 




Ae 




0.034597 




0.087264 


At 






0.166166 


0.104669 


As 


0.070638 


0.247661 






Ag 


0.146412 








Ai4 






0.104300 




Al6 








0.224968 


Al9 






0.114122 




Thresholds 


<5* 


0.480904 


0.481524 


0.491407 


0.491740 



Aa' x^-^ + \l x^-^ with L > 7 and 4 < if < L - 1. For 
any choice of R, dc, K and L, Copt will denote the minimum e 
such that a corresponding p-positive (1 — e) capacity achieving 
distribution of rate R exists, for the given parameters. 
For such A(-), the coefficients of p (•) are: 

Po = (4-2)A2 - (24-2)A3 
Pi = (24-3)A3 

p K-:i = -[{K - l)dc - [K - 1)]\k 

PK-2 = [{K -l)dc-K]XK 

p L-3 = -[{L - l)dc - {L - 1)]Xl 

PL-2 = [{L-l)dc-L]XL. 

According to Corollary |2] necessary (but not sufficient) 
condition for the degree distribution to be p-positive is 
that the number A of sign changes in the sequence 
Po,Pi,Pk-3,Pk-2,Pl-3,Pl-2 is even. In the specific case 
under analysis, the following inequalities are always true: 
PL-2 > 0, pL_3 < 0, pK-2 > 0, pK-3 < 0, pi > 0. Then, 
if and only if p o > 0, the condition that A is even is satisfied 
(specifically, it results A — 4). 

Suppose to fix the code rate R, the check nodes degree 
dc and the active variable degrees (K and L). From the 

constraints X]i^2 ~ ^ i? = 1 — c?c/(X]i^2 '^i/*)' ^^'^ 
from ( fTTT i. expressions of A3 and A^ as functions of e and 
Al, namely A3 — A3(e, A^) and Xk — XK{e,XL), can be 
found. Then, the necessary condition for the distribution to be 
p-positive is given by the following set of inequalities: 

0<A2(e)<l 0<A3(e,AL)<l 
< AA(e,AL) < 1 0<Al<1 
Po(e, Al) > O. 

The solution of this set of inequalities identifies a region on 
the plane e — A^. This region, denoted by A4, will be called 
the permitted region (since no p-positive distribution can exist 
outside this set). Some examples of permitted regions, for 
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Fig. 1. Permitted region for some values of [dc, K, L). 

different values of dc, K and L, all corresponding to R = 1/2, 
are depicted in Fig. [T] If (e, Al) G A4, then the polynomial 
p (•) could have in principle 0, 2 or 4 real roots between 
and 1 (A = 4). 

Let e^jj be the minimum value of e allowed for points 
within A^. If, for e = Cmin, at least a Al exists such that 
(emin, Al) G A^, and such that p{-) is not negative between 
and 1, then it follows eopt = Cmin- On the contrary, if p{x) 
is negative for some x between and 1, for e = emin and for 
any admitted A^, it results eopt > Cmin- This is always the case 
when eniin = 0. An approach to evaluate Copt in this case is 
described next. 

The proposed approach is based on this observation: If 
(copti Ai) is a solution of the optimization problem (for some 
then at least one x G (0, 1) must exist such that 

p(x;eopt, Al) = (12) 
p'(x;eopt, Al) = 0, (13) 

where the dependence of the coefficients of p (■) on e and Al 
have been explicitly indicated. Then, the search for the optimal 
distribution can be restricted to the points (e, Al) G A4 for 
which p (x; e, Al) = and p'(x; e, Al) — 0. 

From these relationships, it is possible to obtain e and Al 
as rational functions of a;, namely e — e{x) and Al = AL(a;). 
Then, the technique consists in plotting the trajectory of the 
point {e{x), Xl{x)) on the plane e — Al, as well as the 
permitted region Ai. From this diagram it is usually possible 
to univocally determine eopt, as shown next with an example. 

For instance, suppose R — 1/2, dc ^ 7, K = 5 and 
L = 17. The trajectory of (e(x), Xl{x)) is depicted in Fig. |2] 
together with a detail of the permitted region. The trajectory 
of (e(x), AL(a;)) crosses the permitted region for values of x 
between xi = 0.527434, and X2 = 0.735514. The minimum-e 
point on the trajectory segment xi X2 is the intersection 
point between the trajectory and the upper boundary of M.. For 
this point, it results e ~ e{xi) ~ 0.032242. The corresponding 
polynomial p{-), whose graph is depicted in Fig. [3] does 
not assume negative values between and 1. Hence, the 
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Fig. 2. Diagram siiowing tlie permitted region and tlie trajectory of tlie point 
(e(s), \n{x)), for R = 1/2, = 7, K = 5 and L = 17. 



P(x) 




0.2 0.4 0.6 0.8 1 

X 



Fig. 3. Plot of p( ) for the optimal R = 1/2 p-positive distribution 
corresponding \o dc = 7 , K = b and L = 17. 

TABLE II 

Optimal R = 1/2 p-positive and not p-positive distributions for 
dc = {6,7} WITH \(x) = A2 X + A3 + A5 x"* + Al x^-^. 





dc = 6, 


L = 10 


dc = 7, L = 15 




p-pos. 


DE 


p-pos. 


DE 


A2 


0.418913 


0.415774 


0.341501 


0.339505 


A3 


0.167565 


0.180916 


0.142292 


0.140214 


As 


0.266696 


0.248100 


0.248395 


0.259036 


Al 


0.146826 


0.155210 


0.267812 


0.261244 


Thresholds 


5* 


0.477426 


0.480325 


0.488041 


0.490947 



minimum-e point on the segment xi X2 corresponds to a p- 
positive distribution: Then, eopt — 0.032242. The threshold of 
the corresponding degree distribution is S* — [A'(0)p'(l)] = 
(l-eopt)(l-i?) = 0.483879, the value of A17 can be obtained 
from the function A17 = Xnix) and the values of A3 and 
A5 can be obtained from the functions A3 = A3(e, A17) and 
A5 = A5(e, A17). 

By adopting a similar approach, we fixed R — 1/2, K = 5, 



and looked for the best p-positive distribution for dc = 6 and 
dc — 7. The optimal distributions were obtained for L = 10 
and L = 15, respectively. They are shown in Table HH as well 
as the optimal check regular distributions obtained by running 
the DE tool, under the same set of constraints (R, dc and 
active variable degrees), but without the p-positive constraint. 
The p-positive thresholds are only slightly lower than the 
not p-positive counterparts, and the degree distributions are 
quite similar. Again, the p-positive distributions have the 
advantage represented by the closed-form threshold expression 
[A'(0)p'(l)]-i. 

V. The Binomial Distributions are ^-positive 

In [7], it was shown that capacity achieving sequences of 
degree distributions for the BEC can be constructed according 
to the binomial degree distribution. This is a check regular 
distribution, whose variable distribution is in the form \{x) — 



"G-i)(-i)' 



, where (^) = a(a-l) . . . {a-N+ 



l)/(iV!), and a = (dc— 1)^^. Recognizing the binomial degree 
distributions as part of the p-positive family, the following 
theorem states that this family can achieve the BEC capacity. 

Theorem 3: Any binomial degree distribution is p-positive. 
Proof: Recall that, for i = 2, . . . , L - 1, the (i - 2)-th 
coefficient of p (•) is pi-2 — oJiXi — {oji+i + 1) A^+i. For the 
binomial distribution, it is 



w.A, = [(4-l)(*-l)-l] 
Furthermore, it is 



{uji+i + 1) A,,+i = {dc - 1) i 



a 



a 



{dc - 1) 



(*-i-«)«G"i)(-ir 



[(4-1) {i-l)-a{dc~l)] 



= [{dc-m-i)-i] 



a-L(2)(-l)^+i 

G"i)(-i)^ 



a-L{f){-l) 



where the second equality is due to the fact that (") — 
ii-i) ""i"'"^ ™d the last equahty to the fact that a ^ {dc - 
Thus, for i = 0, . . . , L — 1, LUiXi — (oj^+i + 1) A^+i, 
i.e. Pi-2 = 0. It follows p {x) = ujl Xlx^^^, which is always 
positive for a; e (0, 1]. □ 

VI. Conclusion 

In this paper, a special family of LDPC degree distributions 
has been presented. The main feature of these distributions is 
the possibility to express in closed form their decoding thresh- 
old on the BEC, under iterative decoding. More specifically, 
the threshold admits the simple closed form [A'(0)p'(l)]^^. 
This family is a subset of the class of check-regular distri- 
butions, and the distributions within this family have been 
called p-positive distributions. A simple necessary condition 



for a check regular distribution to belong to this class has 
been obtained, by invoking some known results about the real 
roots of polynomials. 

Three additional properties of the proposed distributions 
have been highlighted. The first one is their very good 
threshold under some set of constraints. This property is not 
general, depending on the specific imposed set of constraints. 
The second one is the possibility, for particular structures of 
the variable distribution, to optimize the distribution without 
necessarily using the numerical optimization tools which are 
usually exploited. The third one is the possibility to achieve 
the BEC capacity with sequences belonging to the proposed 
family. 
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